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We analyze a decomposition of the Coulomb electron-electron interaction into a long-range and 
a short-range part in the framework of density functional theory, deriving some scaling relations 
and the corresponding virial theorem. We study the behavior of the local density approximation in 
the high-density limit for the long-range and the short-range functionals by carrying out a detailed 
analysis of the correlation energy of a uniform electron gas interacting via a long-range only electron- 
electron repulsion. Possible definitions of exchange and correlation energy densities are discussed 
and clarified with some examples. 



I. INTRODUCTION 

In recent years, there has been a growing interest in ap- 
proaches of density functional theory (DFT) based on 
a long-range/short-range decomposition of the Coulomb 
electron-electron interaction (see, e.g., Refs. 0,0 0,00)- 
The idea is to use different, appropriate approximations 
for the long-range and the short-range contributions to 
the usual density functionals of the Kohn-Sham (KS) 
scheme 0. 

In this paper, we consider one possible long- 
range/short-range decomposition of the Coulombic den- 
sity functionals. To gain insight into the introduced long- 
range and short-range density functionals, we examine 
some scaling relations, the virial theorem and energy den- 
sities. The local density approximation (LDA) appropri- 
ately defined for these functionals is also analyzed. 

The long-range and short-range density functionals 
considered in this work are defined as follows (see also 
Refs. H El IH 13 El The starting point is the 
decomposition of the electron-electron Coulomb interac- 
tion w ee (r) = 1/r as 



w ee (r) 



(r), 



(1) 



with the long-range interaction w^{r) = erf (/xr)/r 
and the complement short-range part w^^ir) — 1/r — 
erf (fir)/r. The parameter /x controls the range of the de- 
composition. For fi — the long-range interaction van- 
ishes, w^~ Q {r) = 0j while for /x — > oo it reduces to the 
Coulomb interaction, mj,'/^ 00 (r) = 1/r. The Coulombic 
universal density functional F[n] is consequently decom- 
posed as 



F[n] =F lr ^[ 



F s 



(2) 



where F h [n] is a long-range universal density functional 
and F sr,A1 [n] is its short-range complement. The long- 
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range functional is defined in the constrained-search for- 
mulation 01 by 



F h 'V[n} = 



min^lf + W / ee M l*) 



(3) 



where f is the kinetic energy operator and = 
(1/2) JJ n 2 (ri, r2)Wg r e ' /i (ri 2 )driffr 2 the long-range inter- 
action operator, expressed with the pair-density oper- 
ator ^(ri,^)- The minimizing (multi-dctcrminantal) 
wave function in Eq. (0 is denoted by \E ,lr,M [n] . In- 
troducing the non-interacting kinetic energy functional, 
T s [n] = ($[n]\f\$[ri\) , where <E>[n] is the KS determinant, 
F lr,/i [n] is written as 



F h ^[n] = T s [n] + E^M 



(4) 



where the long-range Hartree-exchange-correlation func- 
tional, E^ c [n] = E^[n] +E^[n] +E l ^[n], is the sum 
of the long-range Hartree functional 

E n>] =\jj n{v l )n{T 2 )w^t(r l2 )dT l dv 2 , (5) 

the long-range exchange functional 

££■>] = mn]\Wlrmn}) ££•>], (6) 
and the long-range correlation functional 

E x ^[n) = (* h '^[n]|f + W l e l^ h ^[n}) 

-($[n]\f + W%>>\$[n])- (7) 

The short-range functional, F sl '^[n] = E^f c [n] = 
E^[n]+E^[n]+Ef^[n], is the sum of the short-range 
Hartree, exchange and correlation functionals, defined 
by complementarity to the Coulombic Hartree, exchange 
and correlation functionals, £h[w], ^[n] and -E c [rt], 



=E H [n]-E^\n] 



E s x ^{n}=E x {n}-E l ^{n], 



(8) 



(9) 
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£* r <>] = E c [n]-E^[n\. (10) 

The LDA approximation to the long-range exchange- 
correlation functional E*^ DA [n] is constructed from the 

exchange-correlation energy per electron e^tinif °f a un i~ 
form electron gas interacting with potential W"^ 

E l :> c % A [n] = J n{v) e l ^ nnii (n(v)) dr. (11) 

Similarly, the short-range LDA exchange-correlation 
functional B^ldaM ^ s defined with the complementary 

quantity e^ nif = e^cunif - ^cfunif ( see Refs.HIH). 

The paper is organized as follows. In Sec. [HI we exam- 
ine some scaling properties of the long-range and short- 
range functionals, and we discuss the high-density limit 
of the correlation functionals in LDA. In Sec. IIIII we de- 
rive the virial theorem satisfied by the long-range and 
short-range functionals. In Sec. IIVI we examine long- 
range and short-range energy densities. Sec. M contains 
our conclusions. 

All the relations derived in this work are more gener- 
ally true for an interaction of the form vr^(f) ~ (3(/j,r)/r 
where /? is a continuous and differentiable function satis- 
fying f3{x — ► oo) = 1. 

Atomic units (a.u.) are used throughout this work. 

II. SCALING RELATIONS 

In this section, we generalize some usual scal ing: rela- 
tions of the Coulombic density functionals 0, 0] (see 
also Ref. Ils|) to the case of the long-range and short- 
range density functionals. The scaled wave function of 
a iV-electron system corresponding to a uniform scaling 
of the electron coordinate s by the scale factor 7 > is 
defined by (see, e.g., Refs. El E |U H2) 

^(n,...,^) = 7 3A, / 2 *( 7 r 1 ,..., 7 r Ar ). (12) 

The wave function yields the scaled density n 7 

n 7 (r) = 7 3 n( 7 r). (13) 

A. Scaling relation for the Hartree and exchange 
functionals 

It can be easily verified from Eq. JSJ that the long- 
range Hartree functional satisfies the following scaling 
relation 

E^[n 1 ]= 1 E^[n]. (14) 

The KS determinant associated to the scaled density n 7 is 
identical to the KS determinant associated to the density 
n after uniform scaling of the coordinates 

$[n 7 ] = $ 7 [n], (15) 



and thus ($[n 7 ]|W^T|$[n 7 ]) = ($[ n ]|W^|$[n]), lead- 
ing to the same scaling relation for the long-range ex- 
change functional 

E^[n 1 ]= 1 E^[n]. (16) 

The short-range Hartree and exchange functionals satisfy 
the same scaling relations 

E^[n J ]= 1 E^[n}, (17) 



E s ^[n 1 ]= 1 E^[n\. (18) 

B. Scaling relation for the correlation functionals 

Let's introduce the following universal functional 

F ll '"' l [n] = min(¥|f + JW£"|*}, (19) 

and the minimizing wave function is denoted by 
ijjlr^Ajj^ ijv^g sca i ec j wave function ^'^[n] gives the 

density n 7 and minimizes ($\{l/j 2 )f + (A/7)W'ee'* 7 l*) 
since (*i;^' A |(l/7 2 )f + {X h)W^ \^^ x ) 
(^'^if+XWll'^ 11 '^). Consequently, doing the sub- 
stitution A — > A/7, the wave function ^'^'^^[n] gives 
the density n 7 and minimizes (l/ 7 2 )(*|r + AW^ T |#). 
Therefore, we conclude that 

* lr ^- A [n 7 ] = 4^' A /7[ n ]. (20) 

Applying this relation with A — 1, we find the following 
scaling relation for the long-range universal functional 

F h ^[n 7 ] = 1 2 F h ^ 1 ^[n], (21) 

and consequently for the long-range correlation func- 
tional 

£^> 7 ] =7 2 £^< 1 />]. (22) 

where E x *^ l [n] = (* lr ^' A [n]|f + lW^\^ r ^ x [n\) - 
{^[n]\f + lW^ f ; tl \^[n\) . The short-range correlation func- 
tional satisfies the same scaling relation 

Ef'^[ nj ] =7 2 ^f-' i ' 1/7 [n], (23) 

where E s c r ^^[n] = E^[n] - £' r ^ A [n] and £ A [n] is the 
Coulombic correlation functional along the linear adia- 
batic connection. To our knowledge, Eq. I|20|) is new 
while Eq. H21|l has already been mentioned by Yang [2^] • 

C. High-density limit of the correlation functionals 

We study now the long-range and short-range correla- 
tion functionals under uniform scaling to the high-density 



3 



limit (7 — » 00). If the ground-state of the KS system 
is non-degenerate, £ ; ^ r ' M7 [n 7 ] goes to a constant when 

7 — > 00 



lim £' r '^ 7 [nJ = S'^'^Hnl 



(24) 



where 2?, 



lr,M,(2)r 



= (l/2)(d 2 E l *>»> l [n]/dl 2 )i= is the 
second-order correlation energy in the Gorling-Levy per- 
turbation theory [23J, just as in the Coulombic case 0, 
|2~H |. We have a similar behavior for Ef-^lrij] 

lim Ef'^[nJ = E^'^ln], (25) 

7 — ► oo 

with B?'"' (a) [n] = (l/2)(9 2 Sf ^' z [n]/9Z 2 ) z= o. 

It is interesting to study how the long-range and short- 
range LDA functionals behave in the high-density limit, 
and to understand whether they can overcome some 
of the well-known problems of standard LDA in this 
regime (see, e.g., [25|)- For the uniform electron gas of 
density n, the scaling transformation of Eq. i|13[) is sim- 
ply written as r s — > r s /7, where r s — (47rn/3) — - 1 / 3 is the 
Wigner-Seitz radius. The correlation energy per particle 
of the Coulombic electron gas diverges in the high-density 
limit logarithmically .26] 



-A In 7 



(26) 



where A = (1 - ln2)/7r 2 . The Coulombic LDA func- 
tional thus cannot recover the correct high-density scal- 
ing for systems with a non-degenerate ground state |25j . 
The origin of the divergence of Eq. (|2*S|l has been widely 
analyzed in terms of the electron-gas pair density, both 
in real and reciprocal space (see, e.g., [23, IH, HH), 
and, more recently, in terms of the momentum distri- 
bution |3Cj . 

The investigation of the 7 — > 00 performances of the 
long- and short-range LDA functionals implies the study 
of the high-density behavior of the correlation energ y of 
an electron gas interacting with b)^. Following Refs.l27l 
l28l we consider for this "long-range" electron gas the 
Fourier transform of the pair density, the static structure 
factor S '^(k, r s ). The correlation energy per particle, 



Jr,fj, 
"c,unif ' 



is given in terms of this function as 



"c,unif 



dX 



dk 



w^(k)S^' x (k,r s 



(27) 

where S^' X (k,r s ) = S h ^ x (k,r s )- S KS {k,r s ) is the cor- 
relation part of the static structure factor for the system 
of density parameter r s interacting with Xw*^, and 



(28) 



is the Fourier transform of the long-range interaction 
w e^( r ) = erf (/xr)/r. The results derived below hold 
more generally for any long-range interaction whose 
Fourier transform is of the form l|28|) . j£-r](k/fi), with 



rj(x — ► 0) = 1, i.e., for any long-range interaction whose 
real space form is (3{nr)/r with (3(x — > 00) = 1, as stated 
at the end of the Introduction. 

For the Coulomb gas, the random-phase approxima- 
tion (RPA) provides an expression for S x which is ex- 
act for small (<C 1) values of the scaled variable q = 
k/2kp = ar s k/2 (with a 3 = 4/97r) and which gives 
the exact correlation energy when r s — > (see, e.g., 

Refs.lHliJ. The ar S uments for the validity of RPA in 
the high-density limit [3l| can be extended to an inter- 
action of the kind l|28|) as long as, when r s — > 0, w l ^(k) 
diverges for small k as fc~ 2 : in this case, in fact, the 
perturbation series expansion for the correlation energy 
contains as leading term an infinite number of divergent 
direct diagrams whose resummation yields the RPA ex- 
pression for the correlation energy, as in the Coulomb 
gas. The RPA S l ^' x reads 



S. 



\TyfJ,, X 

,RPA 



(<7,?' s ) 



6q 



-v 



2q 



q 2 - \ar s x (q, m) 77 (^§j^) 

(29) 



where xo(q,iu) = (ir /kF)xo(q,iu,r s ) is a reduced non- 
interacting response function, expressed in terms of the 
reduced imaginary frequency u = —iu>/(2qk 2 ? ), that does 
not contain any explicit dependence on r s . In the case 
of the Coulombic gas (which corresponds to fi = 00, that 
is, 77 = 1), if one simply lets r s go to zero, the r.h.s. of 
Eq. H29fl gives a static structure factor which behaves like 
1/q for small q, thus yielding the logarithmic divergence 
of the correlation energy. This divergence comes from 
the combination of the two limits, r s — > and q — ► 0, 
and requires an appropriate treatment, for which it is 
convenient to divide the integral of Eq. (|27|) in two parts, 
Jq 1 dq + J"^° dq, with q\ <C 1. The second part of this in- 
tegral is finite when r s — > 0, while the first part yields the 
logarithmic divergence and it is better studied by switch- 
ing to the scaled variable y = q/y/r^ H3, 0, |3(|, which 
gives an exact scaling for S , ' r ' A1 ' A when q <C 1 and r s — > 0, 



-6 v /r7at/7?(- 



^■"■ A (?«l,r,->0) = 
' \R{u) 2 du 



(30) 



where R(u) = (uarctan^ — l)/ir. Equation H3()|l shows 

that if we rescale //, fj, — ► fi/r s , the factors v( ^Tu^~ ) be- 
come equal to 1 when r s — 0, so that the small-g part of 
g^,n/r a ,\ sca j eg exac t}y ^ ne sam e limit of the Coulom- 
bic gas, i.e., 

5^/ r " A (g<l J r.) ~ -G^F s af(y,X), (31) 

r s — >0 

where the function f(y, A) jlSE^I behaves as y for small y 
and as 1 jy for large y, and is reported in Fig.^in the case 
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FIG. 1: The function f(y, A) of Eq. ijSTjl at A = 1. 



A = 1. Then everything goes as in the Coulomb gas |28j : 
in the small-g part of the integral in Eq. 1)27(1 we can 
replace S^ ,,1 ' r " > ' with Eq. I(31|l . obtaining an expression 

of the kind Jq 1 dyf(y, A). When r s — ► 0, even if qi is 
small, the upper limit of this integral is large, and the l/y 
behavior of f(y ^> 1) causes the logarithmic divergence 
of the correlation energy per particle as 7 — > 00 



Jr,A»7 
'c.unif 



( r s/l) 



-A In 7 + 



(32) 



with exactly the same A of Eq. (|26|l . 

The short-range correlation energy per particle of the 
uniform electron gas is just the difference between the 
correlation energies of the Coulomb gas and of the "long- 
range" gas, so that an exact cancellation of the logarith- 
mic term occurs if we rescale [A, 



lim S(^/7) = h(fj,,r a ). 



(33) 



The function h(/i,r s ) remains finite as long as fir s is 
greater than zero. This means that one can improve the 
LDA performances in the high-density limit by rescaling 
/1 either locally, fj, oc l/r s (r), or globally by choosing a \x 
bigger than l/r™ m , where r™ m is the minimum value of 
r s (r) in the given system. In other words, the rescaled 
short-range functional allows to perform a cut-off of the 
correlation transferred from an electron gas, i.e., we can 
eliminate the long-range correlations that occur in an 
electron gas but that do not occur in a confined system. 



condition implies that the derivative with respect to 7 
vanishes at 7 = 1 



[^[n lh ]\f + Wy\^[n lh ]) 



0. (34) 



7=1 



The kinetic and electron-electron interaction energies 
have simple uniform coordinate scalings 



-( 7 2 T lr <> 1/7 ] + lWlr h [ny y })) ^ = 0, (35) 



where T h ^[n] = (* lr ^[n]|f |* lr >"[n]) and W^[n] = 
(^ lr 'M[ n ]|^lr>M|$lr,M[ n ]^ Performing the derivative with 
respect to 7 leads to 



2T b '>]+W£"[n]-/x 



dW£"[n] 



(36) 



7=1 



Using the virial relation for non-interacting kinetic func- 
tional [3 



(37) 



2T s [n] = ( / • 



7=1 



Eq. i|36|) simplifies to 
T^[n]+E^ c [n]- f 



= - I n(r)r.V ^ H C c [ n] dr. 



— E lt ^\n ' 



5n(r) 



7=1 



(38) 



where T^[n] = T lr ^ [n] - T s [n] and the last equality has 
been obtained through a integration by parts. 

It can be verified that the long-range Hartree, exchange 
and correlation functionals obey separate virial relations, 
just as for the Coulombic case, 



BE, 



lr.jtr 



dfi 



III. VIRIAL THEOREM 



dfj, 



n(r)r.V * . J rfr, (40) 
on(r) 



In this section, we generalize the virial theorem of the 
Coulombic density functionals 0] (see also Refs. |2£j, H2) 
to the case of the long-range and short-range density 
functionals. 

The wave functions of the type \E , ^'^[n 1 / 7 ] give the 
density n, independently of 7, and can therefore be used 
as trial wave functions in the variational definition of 
F h ^[n] [Eq. ©]. As *! y r '' i [n 1/7 ] reduces to the mini- 
mizing wave function ^ lr ' M [n] at 7 = 1, the stationarity 



lrur 1 Irur i dE c^ N f , s SE^Wl] , 

T c i n }+ E c M-f 1 gj~ = - / »(r)r*-^dr. 

(41) 

The same relations are also valid for the short-range 
Hartree, exchange and correlation functionals 
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on J Sn{r) 

(44) 

where T c sr ^[n] = T c [n] - T^[n] and T c [n] is the usual 
Coulombic correlation kinetic functional. 

For the corresponding virial relations in the uniform 
electron gas, see Ref. I^J. 

IV. ENERGY DENSITIES 

In this section, we examine long-range and short-range 
energy densities defined from pair densities or from the 
virial theorem. Energy densities are always useful to an- 
alyze approximations or derive new approximations (see, 

e.g. , Refs. m m m s m m . 

A. Energy densities from pair densities 

In Ref. |4lj, energy densities for modified interactions 
defined from pair densities have been discussed. We re- 
call that an energy density associated to the long-range 
exchange energy E^ can be defined by 

<pd(r) = \J n 2 , x (r,r 12 Wjf(r 12 )dr 12 , (45) 

where n 2 ^ x (r,ri 2 ) is the KS exchange pair density. Like- 
wise, the corresponding short-range energy density writes 

Cp» = \J "2,x(r,r 13 )<'"(r 12 )dr 13 . (46) 

The definition of correlation energy densities from pair 
densities requires an integration over an adiabatic con- 
nection. For instance, an energy density associated to 
the long-range correlation functional E 1 *'^ can be writ- 
ten as 

e c/ P d( r ) = 2J V n 2,c( r . r i2) ^ dr 12 , (47) 

where n 2 r 'c(r, ri 2 ) is the correlation pair density for the 
long-range interaction w^(ri 2 ). The corresponding 
short-range correlation energy density is 

e c >P d( r ) = ^ J * y n 2,£( r > r «) ^ dri2- (48) 

These energy densities involves two-electron quantities 
that can complicate their evaluation. 



B. Energy densities from the virial theorem 

Long-range and short-range energy densities can be de- 
fined from the virial theorem, just as for the Coulom- 
bic case (see, e.g., Refs. [3 Ha)- The virial relation 
of Eq. l(4Tj|l leads indeed to the following long-range ex- 
change energy density 

e^ ir (r) = - / x| °pn(r)r.V^(r), (49) 

where v^(v) = 8E^^[n\/8n(v). Likewise, Eq. leads 
to the following short-range exchange energy density 

e s x %(r) = - M y°°pn(r)r. Vi£>«(r), (50) 

where v x T,tJ/ (r) = SEf'^ln]/ '5n(r). 

The virial relation of Eq. I|4I|) generalized to 
the linear adiabatic connection, T^^[n] + E^'^[n] — 
[idE^inydfM = - fn(r)r.V5E l *^{n]/5n(r)dr where 
T^[n] is given by T^^ l [n) = E^' l [n] - I dE l J'^' l [n]/dl 
(see Refs. 13^, E3) enables to define a long-range correla- 
tion energy density 

/°° r]\ 
-n(r)r-V^(r), (51) 

where v~'^ l ' l (r) = 5E l * ,fJ - l ' l [n]/5n(r). Likewise, the virial 
relation of Eq. I|44|l leads to the short-range correlation 
energy density 

/OO 7 \ 
^n(r)r.V«f^'(r), (52) 

where vf>^(r) = 5Ef^ l [n]/8n(v). 

These energy densities have the advantage of involving 
only one-electron quantities. 

C. Results on the He atom 

As a simple illustration, we have calculated for the 
He atom the long-range and short-range energy densities 

e x,pd( r )> e x,pd( r )> e L r ,vh-M and C'vh'M far /* = 1 a.u., 
as well as the Coulombic energies densities e x , p< j(r) = 

e x,pd( r ) + e x?pd( r ) and e *,virM = e^v h .(r) + e^ r (r). 
For a two-electron system, the exchange pair density and 
exchange potentials are directly deducible from the den- 
sity. Using an accurate density, accurate exchange energy 
densities are thus easily obtained. 

Fig.[2]compares the accurate and LDA radial exchange 
energy densities as a function of the distance to the nu- 
cleus r. One sees that the energy densities defined from 
the exchange pair density and from the virial theorem 
are qualitatively similar. At small r (r < 0.5 a.u.), the 
LDA slightly overestimates the accurate Coulombic en- 
ergy densities. At large r (r > 0.5 a.u.), the LDA im- 
portantly underestimates the accurate Coulombic energy 
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4/rr^e x 
(a.u.) 




47Tr ej 
(a.u.) 



which is represented for the He atom in Fig. [3] with re- 
spect to the interelectronic distance r± 2 for two positions 
of the reference electron r = and r = 1 a.u.. For r = 0, 
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FIG. 3: Accurate (solid curves) and LDA (dashed curves) 
spherical-average exchange hole n| ph ' avr '(r, rw) [Eq. 1531 1 for 
the He atom with r — and r = 1 a.u.. 



FIG. 2: Accurate (solid curves) and LDA (dashed curves) 
radial exchange energy densities denned from pair densi- 
ties (thick curves) and from the virial theorem (thin curves) 
for the He atom: (a) Coulombic radial energy densities 
4nr 2 e x>p a(r) and 47rr 2 e a; , v i r (r) (b) long-range radial energy 
densities 47rr 2 e^ d (r) [Eq. ggj] and 4nr 2 e^ ir (r) [Eq. 
(c) short-range radial energy densities 47rr 2 e^. r ' I ^ d (r) [Eq. 1461 1 
and 47rr 2 e= r ; v A ; r (r) [Eq. jgQfl. with ^ = 1 a.u.. 



both the accurate and LDA exchange holes are centered 
at ri2 = 0, making the LDA a reasonable approximation. 
For r = 1 a.u., the accurate hole is centered near ria = 1 
a.u. while the LDA hole is still centered at r\ 2 = 0, 
leading to an important underestimation of the hole. 



V. CONCLUSIONS 



densities. The contribution at large r remains impor- 
tant in the long-range energy densities while it is signifi- 
cantly reduced in the short-range energy densities. Con- 
sequently, the LDA better performs for the short-range 
energy densities. 

In the case of the energy density defined from the ex- 
change pair density, the better performance of the LDA 
at small distance r can be easily explained in term of 
the spherical average of the exchange hole n x (r, V12) — 
n2, x {r 7 r 12 )/n(r) 

< ph - avl -(r,r 12 ) = ±-j n x (r,r 12 )dn ri2 , (53) 



We have analyzed a short-range and long-range decom- 
position of the Coulomb electron-electron interaction and 
we have derived some exact scaling relations for the corre- 
sponding density functionals. The study of the LDA ap- 
proximation has shown that in the high-density limit the 
short-range functional scales to a constant, thus opening 
the possibility of ameliorating the performances of the 
Coulomb LDA functional in this regime. Possible defini- 
tions of energy densities obtained from pair densities and 
from the virial theorem have been presented. Results on 
the He atom suggest that the LDA approximation can 
give accurate short-range exchange energy densities. 
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